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ABSTRACT 



Current stellar population models have arguably the largest uncertainties in the near-IR wavelength range, partly due to a lack of 
large and well calibrated empirical spectral libraries. In this paper we present a project, which aim it is to provide the first library 
of luminosity weighted integrated near-IR spectra of globular clusters to be used to test the current stellar population models and 
serve as calibrators for the future ones. Our pilot study presents spatially integrated A^-band spectra of three old (>10 Gyr) and metal 

poor ([Fe/H] 1.4), and three intermediate age (1-2 Gyr) and more metal rich ([Fe/H] 0.4) globular clusters in the LMC. We 

measured the line strengths of the Nal, Cal and I2 CO(2-0) absorption features. The Nal index decreases with the increasing age 
and decreasing metallicity of the clusters. The Deo index, used to measure the 12 CO (2-0) line strength, is significantly reduced by 
the presence of carbon-rich TP-AGB stars in the globular clusters with age ~ 1 Gyr. This is in contradiction with the predictions of 
the stellar population models of Maraston (2005). We find that this disagreement is due to the different CO absorption strength of 
carbon-rich Milky Way TP-AGB stars used in the models and the LMC carbon stars in our sample. For globular clusters with age 
>2 Gyr we find Deo index measurements consistent with the model predictions. 
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1. Introduction 

Since the 90's, the interpretation of the integrated light of galax- 
ies (in the nearby universe or at high redshift) heavily relies 
on evolutionary popula tion synthesis (EPS) models. Such mod- 
els were pioneered by Tinsley (1980) and the method has been 
seriously extended since then (e.g . [Bruzual & Chariot! 1993; 
Worthev 1994; Vazdek is et alJI 1 996t iF ioc & Rocca-Volmerange 
Il997t iLeitherer et all 1 19991; IMarastonl 120051; ISchiavonl l2007t) 
They are used to determine ages, element abundances, stellar 
masses, stellar mass functions, etc., of those stellar populations 
that are not resolvable into single stars with today's instrumenta- 
tion, i.e. most of the universe outside the Local Group. To build 
such EPS models we use simple stellar populations (SSP). There 
are two essential advantages of focusing on SSPs. First, SSPs 
can be reliably calibrated. They can be compared directly with 
nearby globular cluster (GC) data for which accurate ages and el- 
ement abundances are independently known from studies of the 
resolved stars. This step is crucial to fix the stellar population 
model parameters that are used to describe model input physics 
and which cannot be derived from first principles (e.g., convec- 
tion, mass loss and mixing). Second, SSPs can be used to build 
more complex stellar systems. Systems made up by various stel- 
lar generations can be modelled by convolving SSPs with the 
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adopted star formation h istory (e.g. iKodama & Arimotol[l997l : 
iBruzual & Charloj [2003). Models describing accurately the in- 
tegrated light properties, including medium to high resolution 
spectra and/or line-strength indices, are and will be our main 
tool to investigate and analyse the star-formation history over 
cosmological time-scales. 

This approach has worked well in the optical spectroscopic 
regime and has led to well calibrated models (e.g . iThomas et al.l 
2003; IBruzual & Charioll200l IMarastonl 120051) . With the ap- 
plication of such models to observed spectra we derive reason- 
able estimates of the main stellar population parameters (age, 
chem ical composition and M/L ratio) in the nearby universe 
(e.g.. lKuntschndl2000l: iTrager et al.ll2000t IThomas et alj |2005h 
ICappellari et al.ll2006t ISanchez-Blazquez et al.ll2007l) as well as 
at higher redshifts (e.g.. iBernardi et al.l 120051; IMarastonl [2005; 
ISanchez-Blazquez et al.ll2009l) . Of course, uncertainties remain 
due to the degeneracy of a ge and metallic ity effects in the opti- 
cal wavelength range (e. g.. IWorthevll 1 9941) . The integrated near- 
IR light in stellar populations with ages > 1 Gyr is dominated 
by one stellar component, cool giant stars, whose colour and 
line indices are m ainly driven by one parameter: metallicity 
(Fro g el et al.| [T978). Near-IR colours and indices also have the 
advantage of being more nearly mass-weighte d, i.e. the near - 
IR mass-to-light ratio is closer to one (see e.g.. lWorthevlll994l) . 
So, by combining the optical as well as near-IR information one 
can resolve the currently remaining degeneracies between age 
and chemical composition, present in the models, and hope to 
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gain a better understanding of star-formation histories. However, 
currently available stellar population models have arguably the 
largest uncertainties in the near-IR and thus it is of paramount 
importance to provide high-quality observational data to validate 
and improve the state-of-the-art models. 

Globular clusters in the Local Group are an ideal laboratory 
for this project since ample information from studies of the re- 
solved stars is available. Yet, integrated spectroscopic observa- 
tions of the Galactic GCs in the near-IR are very challenging due 
to their large apparent sizes on the sky. The Large Magellanic 
Cloud (LMC) and its globular cluster system, located about 50 
kpc away, is a much better observational choice. It shows evi- 
dence for a very complex and still ongoing star formation ac- 
tivity. The LMC GCs have an advantage (for the scope of this 
project) with respect to Galactic GCs - they span a larger range 
in ages. Studies of the LMC globular cluster system show one 
old component with age >10 Gyr. After this time there was a 
"dark age" with just one cluster formed before a ne w burst of 
cluste r formation that has started around 3-4 Gyr ago dDa Costal 
1199 lb . A disadvantage is their lower metallicity. 

The goal of this project is to provide an empirical near-IR li- 
brary of spectra for integrated stellar populations with ages > 
1 Gyr, which will be used to verify the predictions of current 
SSP models in the near-IR wavelength range. Here we present 
the results from a pilot study of Zf-band spectra of 6 globular 
clusters in the LMC. The analysis of their J and //-band spec- 
tra will be discussed in a separate paper. This paper is organised 
as follows: in Sect. |2]we give details about the sample selection 
and the observing strategy. Sect. [3]is devoted to the observations 
and data reduction. In Sect. [4] we discuss the cluster membership 
of the stars in our sample, in Sect. [5] we describe the near-IR 
index measurement procedures. In Sect. [6] we make a compari- 
son between the currently available stellar population models in 
the near-IR with our data, discuss the observed disagreements, 
and give potential explanations. Finally, in Sect. [7] we give our 
concluding remarks. 



2. Sample selection and observational strategy 

Due to our interests in the application of SSP models to the in- 
tegrated light of early-type galaxies, we restricted our sample in 
this pilot study to intermediate age (1-2 Gyr) and old (> 10 Gyr) 
cluste rs. The sample was selected from the B ica et all (1 1996, 
1999) catalogues by choosing the SWB class dSearle et alj l980) 
to be V, VI or VII. In this way we ensured that the target sys- 
tems will have ages > 1 Gyr. We further required the clusters to 
be bright (My < -5.8) and reasonably concentrated (effective 
radius < 35")- Where no literature data were available, the con- 
centration was checked by eye on DSS images. Another selec- 
tion criterion was the availability of auxiliary data, because we 
needed detailed information on age and chemical composition. 
All of th e selected clusters have HST/WFPC2 and/or ACS imag- 
ing (e.g. Olse n et al.lll998HMackev & Gilmorel20 03). integrated 
optical spectroscopy, and spectra of indiv i dual giant stars (e.g . 
Olszewski et alJll99U|Beaslev et alj |2002t Uohnson et alT 2006; 
Mucciarelli et a l. 2008). We can also benefit from near-IR stud- 



ies, both imag ing and spectroscopy, of the giant stars in these 
clusters (e.g. lFrogel et al.llT990T: iMucciarelli et al . 2006), a s well 
as of photometry (Persson et al. 1983; Mucciarelli et al. 2006; 
iPessev et all2006l ). Taking into account the above criteria we se- 
lected six clusters as targets for this pilot project, aiming at val- 
idating the strategy for observations and analysis (see Table Q]). 
Three of the clusters are metal poor (mean [Fe/H] ~ -1.4) and 



have ages of more than 10 Gyr. The other three are more metal 
rich (mean [Fe/H] ~ -0.4) and younger, with ages between 1 
and 2 Gyr. In the literature there are different age and metallic- 
ity estimates for the clusters in our sample, depending on the 
methods used. Here we lis ted the ages based on SWB types, 
given in Frog el et alj £1990) and a compilation of metallicities, 
obtained from the literature. A summary of the clusters' proper- 
ties is given in Table [1] 

Integrated spectra of clusters with less than 10 4 Lq are likely 
to be dominated by statistic al fluctuations i n the number of bright 
AGB and RGB stars (e.g. lRenzinilll998l for more details see 
Sect. |4|. These particular phases of the stellar evolution are one 
of the main contributors to the integrated light of an intermedi - 
ate age stellar population in the near- IR (e.g. iMarastor] 12005). 
To sample as much of the cluster light in a reasonable observ- 
ing time, we made a mosaic of 3 x 3 VLT/SINFONI pointings 
(8" x 8" per pointing with a 0f.'25 spatial sampling) centred on 
each cluster. We estimated the total light sampled by the central 
mosaic (24" x 24") for all clusters using the following equation 
and the 20" radius aperture ZT-band photometry: 



L T = BC K ■ 10 



-0.4-(m K -(m-M)-M 



k.O 



-A K ) 



(1) 



where nix is the observed Zf-band integrated magnitude of the 
cluster, ( m — M) - 18.5 is the adopted distanc e modulus to 
the LMC (Ivan den Bergh|[l998l; lAlves etai1l2002r. I Alvesl l2004t 
Borissova et al. 200 4)), A k is the extinction towards each clus- 
ter dZaritskv et alj|1997l). = 3.33'" is the ZT-band absolute 
magnitude of the Sun dCoxl 120001) . The bolometric correction 
BCk is applied in order to obtain the total bolometric luminosity 
Lj of the cluster. It depends on the adopted age and metallic- 
ity of the stellar populations. In our case it is 0.6 for the group 
of the old and metal poor clusters (NGC 1754, NGC2005 and 
NGC 2019) and 0.3 for the more metal rich and intermediate age 
cluste rs NGC 1806, NGC 2162 and NGC 2173 dMarastonll 19981 
120051) . For all but two clusters we found that the sampled lu- 
minosity is > 10 4 Z<5 (see Table [2}. Only the young clusters 
NGC 2162 and NGC 2173 are at the limit or bellow this value. 

In order to maximise the statistical probability of getting the 
majority of the RGB and AGB stars, we observed, in addition 
to the central mosaics, up to 9 of the brightest stars surround- 
ing each cluster and outside of the central mosaic. Their selec- 
tion was based on K - (J - K) colou r-magnitude diagram s from 
the 2MASS Point Source Catalogue dSknitskie et al.ll20 06) of all 
the stars with reliable photometry and located inside the tidal ra- 
dius o f each cluster (r r taken from McLaug hlin & van der M arel 
I2005L see Table|2|. We selected the stars with (J - K) > 0.9 and 
K > 12.5'" as an initial separation criterion from the LMC field 
population. 

The inclusion of these additional bright stars in our sample 
has a twofold purpose. The original idea, as described above, 
was to provide a better sampling of the integrated cluster light by 
in- or excluding these bright stars (after a careful decision pro- 
cess, based on kinematical and chemical composition assump- 
tions). Second, the stars, which turn out not to be members of 
any cluster, are representative for the LMC field star population 
in the vicinity of our GCs. Thus we obtained an independent 
field AGB stars sample for comparison with the globular clus- 
ters. The main properties of these additional stars are listed in 
TablefJ] 

In Fig.[T]we show the 2MASS Zf -band images, obtained from 
the 2MASS Extended Source Catalogue of our globular cluster 
sample. The black cross and box on each cluster image match the 
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Name 


UT date 


V 


(B-V) 


SWB 


Age 


[Fe/H] 


(1) 


(2) 


(3) 


(4) 


(5) 


(6) 


(7) 


NGC 1754 


2006 Nov 18 


11.57 


0.75 


VII 


10 


-1.42°, -1.54" 


NGC 2005 


2006 Nov 29 


11.57 


0.73 


VII 


10 


-1.35", -1.92 6 , -1.80 c , -1.33 rf 


NGC 2019 


2006 Dec 10 


10.86 


0.76 


VII 


10 


-1.23", -1.18*, -1.37 c ,-1.10 rf 


NGC 1806 


2006 Dec 06 


11.10 


0.73 


V 


1.1 


-0.23 6 , -0.71" 


NGC 2162 


2006 Dec 05 


12.70 


0.68 


V 


1.1 


-0.23", -0.46^ 


NGC 2173 


2006 Dec 06 


11.88 


0.82 


V-VI 


2 


-0.24*, -0.42 / , -0.5 F 



Notes: (1) Cluster name, (2) date of observation, (3) integrated V-b and magnitude, (4) (B - V) colour and (5) SWB type taken from (B ica et al] 



Il996ijl999h . (6) Age of the c luster in Gyr, based on the SWB type dFrogel et alll990h. (7) [Fe/Hl derive d using different methods : "lOlsen et al 
( 1998) - slope of the RGB; * Olszewski et al. ( 1991) - low-resolution Call trip let; c | Johnson et al. I (120061) - high-resolution Fel; d loh nson et al 
( 2006) - high-resolution F e II; " iDirsch et alj d2000l) - Stromgren photometry; ? Grocholski et al. (2006) - low-resolution Call triplet; 
£ |Mucciarelli et alJd2008h - high-resolution spectroscopy. 

Table 2. LMC globular clusters structural and photometric properties. 



Name 


O, 


r, 




^20 


^-'mosaic 


R/nr 


K-MQ6 


(J - K) M 06 


(H - K) M o6 


(1) 


(2) 


(3) 


(4) 


(5) 


(6) 


(7) 


(8) 


(9) 


(10) 


NGC 1754 


1172 


142'.'9 


9.01 


10.21 


2.7 x 10 4 


20" 








NGC 2005 


8'.' 65 


98'.'8 


8.92 


9.80 


3.9 x 10 4 


20" 








NGC 2019 


9772 


121'.' 6 


8.31 


9.21 


6.7 x 10 4 


20" 








NGC 1806 






7.57 


9.19 


3.4 x 10 4 


100" 


7.076 


1.055 


0.271 


NGC 2162 


21'.' 37 


197'.' 2 


9.65 


11.18 


0.5 x 10 4 


120" 


9.071 


1.253 


0.372 


NGC 2173 


34'.' 35 


39375 


8.86 


10.69 


0.9 x 10 4 


60" 


9.050 


1.033 


0.297 



(2005), (4) AT-band integrated magnitude within the tidal radius from Pessev et al. (2006). For NGC 1806, where r, is unknown, we list the 
mag nitude within 200" r adius. The same applies for NGC 2173, because aperture equal to the tidal radius does not exist, (5) K-band magnitude 
from Pessev et al. (2006) with aperture radius 20", which matches reasonably well our central SINFONI mosaics, (6) sampled bolometric 
luminosity within the clusters central mosaics in Lq, computed following Equation Q, (7) integration radius, up to which we have add ed stars 
from the additional bright AGB stars sample, (8), (9) and (10) photometry with fixed aperture radius of 90" of Mucciar elli et alJ ((2006). 



centre and the extent of the SINFONI mosaic coverage, respec- 
tively. The red squares mark the additional bright stars, observed 
around the clusters within the region of the 2MASS image. The 
green circle and cross sh ow the centre and 20" radius aperture, 
which IPessev et all d2006[) used to obtain integrated magnitudes 
for each cluster. These authors have used near-IR images from 
the 2MASS Extended Source Catalogue and have performed 
photometry with different aperture siz es after correcting for the 
extinction and LMC field population. iMucciarelli et alj ((2006) 
have used ESO 3.5 m NTT/SOFI images to provide integrated 
near-IR magnitudes for half of the globular clusters in our sam- 
ple after decontaminating from the LMC field population and 
correcting for the extinction and completeness. The large yel- 
low circles show their 90" fixed radius aperture. We used these 
photometric studies to compare our spectroscopy with integrated 
colours and magnitudes. There is a good agreement between 
the centres of our SINFONI observations and the photometry 
studies, however, in the case of the sparsely populated clus- 
ter NGC 2173 the offset is 17" (in all other cases this offset is 
smaller than 5"). 

A summary colour-magnitude diagram for all observed ob- 
jects in our sample, stars and central mosaics, is shown in 
Fig. [2] For the J and Zf-band magnitude of the ce ntral mosaics 
we ad opted the 20" radius aperture photometry of IPessev et alj 
(2006), which, as shown on Fig.Q] matches reasonably well our 
central mosaics. The selected in addition bright stars outside 
the central mosaics are denoted with diamond symbols. Their 
photometry comes from the 2MASS Point Source Catalogue 
an d magnitude s were dereddened following the same method 
as IPessev et al] d2006l) - extinction values were obtained from 
the Magellanic Clouds Photometric Survey (Zarit skv et alj 1997b 



and adopting the extinction law of Bes sell & Brett! (Q988). 
The slanted line in this figure represents the separat i on be - 
tween oxygen- and carbon-rich giant stars of lCioni et al. (2006). 
According to this criterion we have five carbon rich stars in our 
sample. 

3. Observations and data reduction 

3.1. Observations 

The observations of the selected globular clusters and stars were 
obtained in service mode in the period October - December 
2006 (Prog. ID 078.B-0205, PI: Kunts chner). We used the in- 
tegral field unit spe ctrograph SINFONI (Eisenh auer et al.f 2003; 
Bonnet et al. 2004), which is mounted in the Cassegrain fo- 
cus of Unit Telescope 4 (Yepun) on VLT at Paranal La Silla 
Observatory. Its gratings are in the near-IR spectral domain (1 
- 2.5 jt/m). We used the Zf-band grating, which covers the wave- 
length range from 1.95 to 2.45 /mi at a dispersion of 2.45 A/pix. 
The spectral resolution around the centre of the wavelength 
range is R ^ 3500 (6.2 A FWHM), as measured from arc 
lamp frames. The largest single pointing with SINFONI covers 
8" x 8", which is too small even for the most compact LMC 
clusters due to their relatively large apparent sizes on the sky. To 
sample at least one effective radius of our targets we decided to 
use a 3 x 3 mosaic of the largest FoV of SINFONI, thus covering 
the central 24" x 24" for each cluster. Given the goal to sample 
the total light, and not to get the best possible spatial resolution, 
all the observations were performed in natural seeing mode, i.e. 
with no adaptive optics correction. The integration time was cho- 
sen based on the requirement to achieve a signal-to-noise ratio of 
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Fig. 1. K-band 2MASS images of our cluster sample. The black boxes and crosses represent our SINFONI 24" x 24" mosaic FoV 
and centres, respectively. The red squares mark the additional bright stars within the region of the 2MASS images that we have 
obse rved around each cluster. The green circles and crosses represent the 20" radius aperture and the centre used in the photometry 
bv lPessev et alJd2006l). The large ye llow circles on the images of the intermediate age clusters (bottom row) illustrate the 90" radius 
aperture of Mucciarelli et al. (2006). The size of each image is 3.'5 x 3.'5. North is up, east - to the left. The very bright star to the 
south-east of NGC 1754 was identified as galactic foreground star and therefore was not included in the analysis. 



at least 50 in the final integrated spectra. The exposure time for 
one pointing of the mosaic was 150 s, divided into three integra- 
tions of 50 s, dithered by 0725 to reject bad pixels. With the short 
integration time we could reliably measure and subtract the very 
bright and variable near-IR night sky using the data reduction 
procedure as described in the next section. 



In order to correct for the effects of the night sky, we ob- 
served empty sky regions very close in time and space to our 
scientific observations, in a SOOOS sequence for each mosaic 
pointing (S - sky integration, O - object). For the additional 
bright stars outside the central mosaic we used the same se- 
quence, but with shorter integration times of 10 s per individual 
integration, leading to a total on source time of 30 s. The sky 
fields for each cluster were located outside its tidal radius (see 
Table |2]l and were checked by eye on 2MASS images to be de- 
void of bright stars. The total execution time for the longest ob- 
serving sequence did not exceed 1.5 hours. This ensured that the 
telluric correction, derived from the telluric stars, observed after 
each cluster and additional bright stars sequence (see Table |4|i 
would be sufficiently accurate. The telluric stars were observed 
at similar airmass as the clusters. 



3.2. Basic data reduction 

An overview of the different data reduction steps and their re- 
sults is shown in Fig. [3] There we show a sequence of a raw 
spectrum, then the sky subtracted spectrum, the telluric correc- 
tion spectrum, and the final, fully reduced s pectrum. More de - 
tails about the data reduction can be found in lLvubenoval (12009). 
Here we briefly discuss the most important steps. 

The basic data reduction was performed with the ESO 
SINFONI Pipeline v. 1.9.2. Calibration products such as distor- 
tion maps, flat fields and bad pixel maps, were obtained with the 
relevant pipeline tasks ("recipes"). To reduce the clusters and ad- 
ditional star data, we divided each observing sequence into clus- 
ter frames (27 object plus 10 sky exposures) and star frames (3 
star exposures plus one sky per star). This was needed due to the 
different integration times for these two sub-sets. We also pre- 
ferred to reduce each mosaic pointing separately and combine 
the nine later. In this way we controlled the quality of each on- 
source frame sky correction and where needed we tuned some 
of the parameters in the pipeline. In summary, we have fed the 
sinf o_rec_j itter recipe with data sets, consisting of one mo- 
saic pointing (3 on source integrations plus 2 bracketing "skies") 
and the needed calibration files. Then the recipe extracts the 
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Name 




R. A. 


Dec. 


K 


(J-K) 




Deo 


Notes 


Cluster 


(1) 




(2) 


(3) 


(4) 


(5) 




(6) 


(7) 


(8) 


2MASS J04540 127-7026341 


04 


54 01.27 


-70 26 34.17 


11.40 


1.13 


1 


,208±0.013 




NGC 1754 


2MASS J04540771-7024398 


04 


54 07.72 


-70 24 39.88 


11.86 


1.26 


1 


,223±0.021 






2MASS J04543522-7027503 


04 


54 35.22 


-70 27 50.36 


12.26 


1.07 


1 


,143±0.015 






2MASS J04542864-7026142 


04 


54 28.64 


-70 26 14.23 


12.62 


0.91 


1 


.112+0.018 






2MASS J04540935-7024566 


04 


54 09.35 


-70 24 56.66 


12.73 


0.96 


1 


.088+0.021 






2MASS J04541 188-7028201 


04 


54 11.88 


-70 28 20.13 


13.04 


0.91 


1 


,177±0.030 






2MASS J04540536-7025202 


04 


54 05.36 


-70 25 20.20 


13.36 


0.93 


1 


,122±0.029 






2MASS J05 30222 1-6946 124 


05 


30 22.21 


-69 46 12.48 


10.65 


1.21 


1 


,258±0.008 




NGC 2005 


2MASS J05300708-6945327 


05 


30 07.08 


-69 45 32.73 


11.66 


1.04 


1 


,213±0.012 






2MASS J05295466-6946014 


05 


29 54.66 


-69 46 01.44 


11.68 


1.12 


1 


,246±0.010 






2MASS J05 300704-694403 1 


05 


30 07.04 


-69 44 03.11 


11.80 


1.03 


1 


,251±0.010 






2MASS J05 300360-69443 1 1 


05 


30 03.60 


-69 44 31.18 


12.02 


0.97 


1 


,216±0.013 






2MASS J05295822-6944445 


05 


29 58.22 


-69 44 44.59 


12.05 


1.22 


1 


,233±0.016 






2MASS J05320670-7010248 


05 


32 06.70 


-70 10 24.84 


10.49 


1.13 


1 


,136±0.014 




NGC 2019 


2MASS J05 3 1 3 862-70 1 0093 


05 


31 38.62 


-70 10 09.35 


10.56 


1.39 


1 


,196±0.017 


C 




2MASSJ053 15232-7010083 


05 


31 52.32 


-70 10 08.39 


10.61 


0.94 


1 


,200±0.016 






2MASSJ05321 152-7010535 


05 


32 11.52 


-70 10 53.52 


10.87 


1.06 


1 


,176±0.018 






2MASS J05321647-7008272 


05 


32 16.47 


-70 08 27.26 


10.95 


1.13 


1 


,183±0.016 






2MASSJ05320418-7008151 


05 


32 04.18 


-70 08 15.15 


11.22 


1.02 


1 


,204±0.024 






2MASS J0502 1232-6759369 


05 


02 12.32 


-67 59 36.92 


10.23 


1.84 


1 


,051±0.013 


+,c 


NGC 1806 


2MASS J05020536-6800266 


05 


02 05.36 


-68 00 26.69 


10.63 


1.61 


1 


,139±0.015 


+,c 




2MASS J05015896-6759387 


05 


01 58.96 


-67 59 38.76 


10.69 


1.76 


1 


.065+0.013 


+,c 




2MASS J0502 1623-6759332 


05 


02 16.23 


-67 59 33.22 


11.02 


1.06 


1 


.211+0.012 


+ 




2MASSJ05021870-6758552 


05 


02 18.70 


-67 58 55.20 


11.32 


1.11 


1 


.267+0.015 


+ 




2MASS J0502 1 846-6759048 


05 


02 18.46 


-67 59 04.82 


11.74 


1.00 


1 


.243+0.018 






2MASS J0502 1 1 2 1 -6759295 


05 


02 11.21 


-67 59 29.54 


11.97 


0.96 


1 


.197+0.026 






2MASS J0502 1 1 37-675 840 1 


05 


02 11.37 


-67 58 40.15 


11.98 


1.02 


1 


.199+0.023 






2MASS J06002748-6342222 


06 


00 27.48 


-63 42 22.29 


9.60 


1.80 


1 


.075+0.012 


+,c 


NGC 2162 


2MASSJ06003156-6342581 


06 


00 31.56 


-63 42 58.14 


11.64 


1.03 


1 


.202+0.011 


+ 




2MASS J060033 1 6-6342 1 3 1 


06 


00 33.16 


-63 42 13.18 


12.24 


0.99 


1 


.197+0.016 


+ 




2MASS J06003 869-634 1393 


06 


00 38.69 


-63 41 39.37 


12.26 


0.99 


1 


.184+0.020 






2MASS J05563892-7258155 


05 


56 38.92 


-72 58 15.52 


11.77 


1.03 


1 


.227+0.012 




NGC 2173 


zMAa!) JU55 1 566 /- /2jo299 


U5 


a zz. £~h 
5 / 56.6/ 


ti on n i 
-11 58 2y.91 


12.0/ 


0.95 


i 
1 


.160+0.014 


+ 




2MASS J05570233-7257449 


05 


57 02.33 


-72 57 44.95 


12.13 


1.04 


1 


.159+0.016 






2MASS J05575784-7257548 


05 


57 57.84 


-72 57 54.89 


12.18 


1.03 


1 


.203+0.016 


+ 




2MASS J05563368-7257402 


05 


56 33.68 


-72 57 40.28 


12.43 


1.00 


1 


.146+0.016 






2MASSJ05581 142-7258328 


05 


58 11.42 


-72 58 32.87 


12.45 


1.04 


1 


.154+0.019 


+ 




2MASS J05583257-7258499 


05 


58 32.57 


-72 58 49.90 


12.48 


0.92 


1 


.136+0.015 






2MASS J05572334-7256006 


05 


57 23.34 


-72 56 00.66 


12.56 


1.01 


1 


.227+0.019 






2MASS J05565761-7254403 


05 


56 57.61 


-72 54 40.38 


12.84 


1.01 


1 


.169+0.020 







Notes: (1) 2MASS catalogue star name, (2) and (3) star coordinates (J2000) given in hours, minutes and sec onds, and degrees, arc min and arcsec, 
(4) extinction corrected /f-band magnitude and (5) (J - K) colour from the 2MASS Point Source Catalogue ( Skrutskie et al. 2006(), (6) Deo index 
value, measured from our spectra, (7) Notes on individual stars: "C" - a carbon-rich stars, "+" - the star was used for the integrated spectrum of 
the cluster, (8) globular cluster, next to which the star was observed. 

raw data, applies distortion, bad pixels and flat-field corrections, 
wavelength calibration, and stores the combined sky-subtracted 
spectra in a 3-dimensional data cube. The same pipeline steps 
were also used to reduce the observations of the additional bright 
stars and telluric stars. 

The main difficulty during the sky correction arises from 
the fact that our observations were sky dominated, combined 
together with the small amount of flux in some of the frames. 
For these cases we found that by appropriately setting the pa- 
rameters that indicate the edges of the object spectrum loca- 
tion (--skycor-llx , --skycor-lly , --skycor-urx and 
--skycor-ury), we could achieve a good sky correction. 

3.3. Telluric corrections 

The next data reduction step was to remove the absorption fea- 
tures originating in the Earth's atmosphere. These features are 



especially deep in the blue part of the /if -band (blue-wards of 
2.1 //m). For this purpose we observed after each science target 
a telluric star, which is of hotter spectral type (usually A - B 
dwarfs, see Table HI). Since these stars are hot stars, we know 
that their continuum in the A'-band is well approximated by the 
Rayleigh- Jeans part of the black body spectrum, associated with 
their effective temperature. They show only one prominent fea- 
ture, the hydrogen Bracketty absorption line at 2.166//m. For 
each telluric star spectrum first we modelled this line with a 
Lorentzian profile, with the help of the IRAF task splot, and 
then subtracted the model from the star's spectrum. Then we di- 
vided the cleaned star spectrum by a black body spectrum with 
the same temperature as the star to remove its continuum shape. 
Thus we obtained a normalised telluric spectrum. The last step 
before applying it to the science spectra was to scale and shift 
in the dispersion direction each telluric spectrum for each data 
cube by a small amount (< 0.5 pix, 1 pix = 2.45 A) to min- 



6 



M. Lyubenova et al.: Integrated A"-band spectra of globular clusters in the LMC 



Table 4. LMC telluric stars observing log. 
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Fig. 2. Colour-magnitude diagram including all the objects 
in our observational sample. The photometry of the six GCs 
(coloured square symbols) comes from the catalogue by 
Pess ev et"ai1 (|2006). Data for the additional bright stars (dia- 
mond symbols) come from the 2MASS Point Source Catalogue 
dSkrutskie et al.l 12006b . Colour coding of the symbols for the 
bright stars matches the cluster in whose vicinity they were ob- 
served. Stars associated with the old clusters are likely not clus- 
ter members as discussed in Sect. 14.11 The slanted line shows 
the separation betw een oxygen-rich (le ftwards) and carbon-rich 
(rightwards) stars of lCionietaTI(l2006l) . 

NGC2173 



6 - 



g 4 



-Q 
< 



Sky subtracted spectrum 



Telluric spectrum 



Fully reduced spectrum 



2.0 2.1 2.2 2.3 

Wavelength [urn] 



2.4 



Fig. 3. Overview of the data reduction steps, applied to obtain 
fully reduced spectra (in this case, the LMC globular cluster 
NGC2173). From top to bottom: (1) Raw spectrum, (2) Sky 
subtracted spectrum after running the SINFONI pipeline. (3) 
Telluric correction, used to remove the telluric absorption fea- 
tures. (4) Fully reduced spectrum, used for analysis. All data re- 
duction steps were performed to the full data cubes and then, 
after the telluric correction, integrated spectra were derived. 



imise the residuals of the telluric lines (for more details about 
this procedure see Silva et al. 2008). After that each individual 
cluster mosaic and star data cube was divided by the so opti- 
mised telluric spectrum. In this way we achieved also a relative 
flux calibration. 

One telluric star, HD 44533, used for the telluric correction 
of NGC 2019 and its surrounding stars, has an unusual shape of 



Name 


Cluster 


K. A. 


Dec. 




(i) 


(2) 


(3) 


(4) 


(5) 


HD 40624 


NGC 1754 


05:54:56.97 


-65:53:00.65 


A0V 


HD43107 


NGC 2005 


06:08:44.26 


-68:50:36.27 


B8V 


HD 44533 


NGC 2019 


06:14:41.94 


-73:37:35.44 


B8V 


HD 42525 


NGC 1806 


06:06:09.38 


-66:02:22.63 


A0V 


HD 45796 


NGC 2162 


06:24:55.79 


-63:49:41.32 


B6V 


HD 46668 


NGC 2173 


06:27:13.64 


-73:13:42.08 


B8V 



Notes: (l)Telluric star name, (2) Cluster for which this star was used, 
(3) and (4) Coordinates of the star, as listed in Simbad (July 2008), (5) 
Spectral type. 



the Brackett y line. It seems to have also some emission together 
with the absorption. To remove it, we interpolated linearly the 
region between 2.1606 and 2.1706 fim. In this region there are 
not many strong telluric lines, but this interpolation will reflect 
in an imperfect correction of the science spectra for this cluster 
at the above wavelengths. However, none of the spectral features 
of interests for this study reside in this wavelength range. 



3.4. Cluster light integration 

As a result from the previous data reduction steps we obtained 
fully calibrated data cubes for each SINFONI pointing, where 
the signatures of the instrument and the night sky are removed 
as well as possible. During the next step we reconstructed the 
full mosaic of each cluster. For example, Fig.|4]shows the recon- 
structed image of the cluster NGC 1754, together with a 2MASS 
A'-band image for comparison. In this image we still see the im- 
prints of the edges of the individual mosaic tiles. On both im- 
ages, the 2MASS and even better on the SINFONI image, we 
see individual stars, which can be extracted and studied sepa- 
rately. 

However, here we are interested in luminosity weighted, in- 
tegrated spectra for each cluster, to compare with stellar popula- 
tion models. To construct one integrated spectrum per cluster we 
first estimated the noise level in each reconstructed image from 
the mosaic data cube. We considered this noise is due to residu- 
als after the sky background correction. Thus we computed the 
median residual sky noise level and its standard deviation, after 
clipping all data points with intensities more than 3<x (assuming 
that these are the pixels, which contain the star light from the 
cluster). After that we selected all spaxels, which have intensity 
more than three times the standard deviation above the median 
residual sky noise level. We summed them up and normalised 
to 1 s exposure time. In some of the spectra (NGC 1754 and 
NGC 2019), we still suffered from sky line residuals, originat- 
ing from the addition of imperfectly sky corrected spaxels (this 
may happen with intrinsically low intensity spaxels). In these 
cases we interpolated the contaminated regions. 

Our observations were carried out in service mode with con- 
straint sets allowing seeing up to 2". For the individual stars this 
led to a failure of the standard pipeline recipe while extracting 
ID spectra. Moreover in some cases in the field-of-view there 
were also other stars. Thus we decided to manually control the 
selection of star light spaxels. We used the same method as for 
the central clusters mosaics to obtain the spectra of the additional 
stars. 
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Fig. 4. Two views of NGC 1754 in the jf-band. The left image 
is fro m the 2MASS Extended Source Catalogue ( Skr utskie et al.l 
120061) . The white square marks the field, which our SINFONI 
mosaic observations cover (24" x 24"). The right image is re- 
constructed from the data cube, built up from the SINFONI ob- 
servations. The very bright star to the south-east of the central 
cluster mosaic seen in the left image was identified as a galactic 
foreground star. 



3.5. Error handling 

The SINFONI pipeline does not provide error estimates, which 
carry information about the error propagation during the differ- 
ent data reduction procedures. Thus we have to use empirical 
ways of computing the errors. For this purpose we derived a 
wavelength dependent signal-to-noise ratio (S/N) for each clus- 
ter integrated spect rum using the empirical method described by 
IStoehr etall d2007l) . We computed the S/N in 200 pixels width 
bins from the science spectra (corresponding to 0.049yum wave- 
length intervals). Then we fitted a linear function to the S/N val- 
ues in the range 2. 1 - 2.4 yum where all of the spectral features of 
interest reside. The error spectrum for each cluster and star was 
derived by dividing the science spectrum by the so prepared S/N 
function. The selection of the bin width was made after experi- 
menting with a few smaller and larger values. In the case of bin 
widths of 50 or 100 pixels (0.01225 or 0.0245 fim, respectively), 
the S/N function becomes very noisy. Choosing wider, 300 or 
400 pixels bins flattens the features of the S/N function. In gen- 
eral the S/N decreases with increasing wavelength. This is due 
to the combination of the spectral energy distribution and the in- 
strument + telescope sensitivity. The mean S/N around 2.3 fim is 
larger than 50 for the integrated spectra over the central mosaics 
of the clusters and in the range 10-40 for the individual stars, 
depending on their magnitude. The S/N estimate for carbon star 
spectra is a lower limit due to the numerous absorption features 
due to carbon based molecules, which are interpreted by the rou- 
tine as noise. The error spectra were used to estimate the errors 
of the index measurements, as shown in the plots in this paper. 

4. Additional bright AGB stars - cluster members or 
not? 

As discussed in Sect.|2]one of the main problems, when one tries 
to compile a representative integrated spectrum of a globular 
cluster, is the stochastic sampling of bright AG B stars for clus- 
ters with modest luminosities. As an example, iRenzin il (119981) 
shows that a stellar cluster with 1O 5 L0, solar metallicity and age 
of 15 Gyr will have about 1200 red giant branch stars (RGB), 30 
early asymptotic giant branch stars (E-AGB), and only 2 ther- 
mally pulsating asymptotic giant branch (TP-AGB) stars. The 
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latter stellar evolutionary phase is particularly important for in- 
termediate age stellar populations of ~ 1 Gyr, s ince up to 80% 
of their total ^-band light originates there (e.g. Maraston 2005, 
and references therein). In general the old clusters (NGC 1754, 
NGC 2005 and NGC 2019) are massive and well concentrated, 
with half-light radii sampled with our SINFONI central mosaics. 
Therefore our observations sample > 50% of the total bolomet- 
ric light coming from these clusters, and in all cases we sample 
at least ~ 3 - 7 x 10 4 L o . 

However, we were not performing that well with the inter- 
mediate age clusters, where for two of them, NGC 2162 and 
NGC 2173, we sample even less than 10 4 Lq (see Table|2]). Our 
observations are intrinsically affected by statistical fluctuations 
in the number of bright stars, due to the modest luminosities of 
the young clusters. These clusters are less massive and less con- 
centrated than other clusters in our sample. In the more massive 
clusters, mass segregation has likely caused a central concen- 
tration of the more massive main sequence progenitors of the 
observed AGB stars than the lower mass background stars that 
have not yet evolved off the main sequence. In the younger, less 
massive clusters, mass segregation has likely been less efficient. 
Hence, the AGB progenitors will be less centrally concentrated 
relative to other unevolved stars. In turn, this implies that AGB 
stars may be found at larger projected radii in lower mass clus- 
ters than in higher mass clusters. 

In order to have integrated spectra as representative of our 
globular clusters, including the most important stellar evolution- 
ary phases, we observed a number of additional bright stars with 
near-IR colours and magnitudes in the range expected for RGB 
and AGB stars, as explained in the previous sections. After prob- 
ing their membership, we added some of them to the central mo- 
saics cluster light to obtain the final integrated spectra, as ex- 
plained bellow. We have to note that we have added only integer 
numbers of bright stars. This is useful when one aims to ob- 
tain a representative spectrum for a given globular cluster, what 
was our goal in this project. In order to achieve an integrated 
spectrum representative of a full stellar population with a given 
age and chemical composition, where stochastics does not play a 
role, one should also consider adding fractional number of bright 
AGB stars. 

The first criterion for cluster membership of the stars is the 
proximity to the cluster centre, but a LMC field star might also 
have a relatively nearby position due to projection effects. The 
second possibility is to explore the radial velocities of the nearby 
stars and the cluster under study. For this purpose we also need 
to know the observed velocity dispersions of the stars in the clus- 
ters. Due to the low velocity resolution of our observations we 
used data from the literature. For the old clusters (NGC 1754, 
N GC 2005 and NGC 2 019) we took these values from the study 
of iDubath et al.l dl997l) . They have measured core velocity dis- 
persions from integrated optical spectra, covering the central 
5" x 5" of each cluster. The values are 7. 8+ 3.0km s -1 for 
NGC 1754, 8.1 ± 1.3 km s" 1 for NGC 2005 and 7.5 ± 1.3 km s" 1 
for NGC 2019. For the intermediate-age clusters we could not 
find similar observed velocity dispersions in the literature, thus 
we used the predicted li ne-of-sight velocity dispersio n at th e 
centre of the cluster by iMcLaughlin & van derMarell (120051) . 
The numbers are 1.1 km s" 1 for NGC 2162 and 2.0km s" 1 for 
NGC 2173. We did not find a similar estimate for NGC 1806, so 
we used a conservative upper limit of 8 km s . 
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4.1. Old clusters 

In the case of the three old clusters, NGC 1754, NGC 2005 and 
NGC2019 we cannot exclude any star around any cluster, be- 
cause their velocities are consistent with cluster membership. 
This is due to the insufficient velocity resolution of our obser- 
vations. However, some of the brightest stars that we have ob- 
served are located closer to the tidal radii of the clusters, than to 
their centres. This is not expected for old clusters, which have 
already undergone a core collapse and have their most mas- 
sive stars concentrated to wards the centre of the cluster. Indeed, 
Mac kev & Gilmorel (120031) classify the old clusters in our sam- 
ple as potential post-core-collapse clusters, due to the very well 
expressed power-law cusps in their centres. Moreover, in Fig. [2] 
the (J - K) colours of the central mosaics are much bluer than 
the colours of the additiona l bright stars around the old clus- 
ters. ISantos & Frogell (1 19971) point out that very young and not 
massive clusters can be systematically bluer than average due to 
stochastic sampling of the IMF. However, the integrated spec- 
tra of the three old globular clusters discussed here are sam- 
pling several times 10 4 L (T) and have ages of ~ 1 Gyr. In this 
regime the simulations of ISantos & Frogell (1 19971) predict much 
less fluctuations in the {J - K) colour than the observed differ- 
ence between the central mosaics and the additional stars. This 
additionally led us to the conclusion that these stars are not likely 
to be cluster members. In the following we considered them as 
members of the LMC field population. The bright star marked 
with a red square just outside the SINFONI field-of-view for 
NGC 1754 in Fig. [T]might be a cluster member, but the S/N of 
its spectrum is too low (~ 10) and adding it to the final integrated 
cluster spectrum would not increase the total S/N. 

4.2. Intermediate age clusters 

This sub-sample includes the clusters NGC 1806, NGC 2162 and 
NGC 2173. The last two of them are poorly populated as seen in 
near-IR light, visible from Fig. Q] so in their cases the poten- 
tial inclusion of additional bright stars to the central mosaic is 
very important for the total cluster light sampling. A detailed 
photometric study o f the RGB and AGB stars in these clus- 
ters is available from lMucciarelli et al.l (2006). Based on near-IR 
colour-magnitude diagrams and after decontamination from the 
field population, they report several stars in these evolutionary 
phases, as well as their luminosity contribution to the total light 
of the clusters. They consider all stars brig hter than K » 12.3 , 
which represents the level of the RGB tip dFerraro et al . 2004), 
and (J - K) between 0.85 and 2.1 to be AGB stars. 

The AGB stars separate in oxygen rich (M-stars) and carbon 
rich (C-stars). An AGB star becomes C-rich, when the amount of 
dredged up carbon into the stellar envelope exceeds the amount 
of oxygen. Then all the oxygen is bound into CO molecules. The 
remaining carbon is used to form CH, CN and C2 molecules. 
This process is more effective in more metal poor ste llar popu- 
lation s, thus they are expected to have more C-stars ( Ma rastonl 
120051) . C-stars can cont ribute up to 60% of the total luminosity 
of metal-poor clusters (jFrogel et alJ[l990h . Another important 
statement that Frog el et ail (11990) make is that C and M type 
stars are found both in clusters and in the LMC field. Thus it is 
very difficult to separate intermediate age globular cluster stars 
from the field population. 

However, the LMC field carbon stars contamination is not 
expected to be significant at the locations of the globular clusters 
in our sample. The upper limits, according to the carbon stars 



frequency maps of IBlanco & McCarthy] ([1983) are of the order 
of 0.05 to 0.7 C-type stars for an area with a radius of 100". 

NGC 1806: This cluster is relatively rich in stars, as seen on 
Fig. Q] However, there is a lack of information about its spec- 
tral properties. So far, studies of the stellar popula tions have 
been made mainly with photometric methods (e.g. [F rogel et al. 
1990; iDirsch et al.ll2000t iMucciarelli et al.ll2006t iMackev et al. 
2008 ), and two brigh t cluster stars have been observed by 
Olsz ewski et al.l (Il99lh to spectroscopically estimate the metal- 
licity of the cluster (see Table [TJ. Moreover, there is neither in- 
formation about its dynamical properties in the literature, nor 
information about its bounds. Thus we had to chose an arbitrary 
limiting radius for the bright stars inclusion. 

Muc ciarelli et all {2006) point out the presence of 75 RGB 
stars and 9 AGB stars, from which 4 are carbon rich stars, in a 
radius of 90" from the centre of the cluster. Indeed, investigating 
the spectra of the resolved stars in our SINFONI central mosaic 
data cube for this cluster, we identified one of the stars as C- 
type. The first three of the brightest additional stars within 90" 
are also C-type. T hus we chose to integrat e all the stars within 
the radius used by Muc ciarelli et al.l (|2006), with the exception 
of the three faintest stars due to their low S/N spectra. 

NGC2162: Inside the ap erture photometry radius of 90", 
that IMucciarelli et al.l d2006) have used, there is a very bright 
carbon star with K = 9.60 m and (J - K) = 1.80 (see Table©. 
It is visible in Fig. Q] at ~50" northwest of the cluster centre. If 
this star is a cluster member, it would be responsible for ~ 60% 
of the total /f-band cluster light and will affect significantly the 
integrated spectral properties of the cluster, thus it is very impor- 
tant to carefully evaluate its cluster membership. The velocity 
of this star is within the errors the same as the velocity, which 
we measured from the central cluster mosaic. However this is 
not a definitive proof of its membership, as discussed above. 
Looking at the surface distribution maps o f C- and M-type gi- 
ants a cross the LMC field in the study of IBlanco & McCarthy! 
(119831) . we see that NGC 2162 is located in a region far away 
from the LMC bar, which has the highest frequency of field 
carbon stars. According to these maps, we can expect to have 
~ 25 C-type stars in an area of 1 deg 2 . This density, scaled to the 
area, covered by a circle with a radius of 90", gives 0.05 carbon 
stars for the 1 .96 x 10~ 3 deg 2 . Moreover, the control field that 
IMucciarelli et al. (2006) used to estimate the LMC field contri- 
bution (shown on their Fig. 4), does not contain any stars on the 
AGB redder than {J - K) = 1.2. Having bright carbon stars is 
not untypical for LMC globular cl usters with the age a nd metal- 
licity of NGC 2162, as shown by iFrogel et all < fl990b . With all 
these facts taken together, we conclude that this very bright and 
red star is most probably a cluster member, although its definite 
membership will be confirmed or rejected only by high resolu- 
tion spectroscopy. In order to obtain the final integrated spectrum 
of NGC 2162, we have also added two more bright stars, located 
within 90" radius. The remaining stars in our observational sam- 
ple are fainter and their spectra have too low quality. Thus they 
would not increase the total S/N. 

NGC 2173: Deciding about the cluster membership of the 
stars observed around this cluster was easier, due to the availabil- 
ity of high resoluti on spectroscopy of RGB stars in NGC 2173 
from the study of Mucciarell i et al.l (120081) . As shown in Fig. [5] 
we have three stars outside of and one within the central 
SINFONI mosaic i n common with their s tudy. Based on high 
resolution spectra IMucciarelli et al.l d2008) find that these four 
stars have very similar radial velocities (rms = 1.2km s -1 ) and 
negligible star-to-star scatter in the [Fe/H]. Thus we safely con- 
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Fig. 5. A summary of the available photometric and spectro- 
scopic data about NGC 2173. The white square marks the extent 
of the SINFONI mosaic. The red squares show the additional 
bright stars we have observed with SINFONI. The cyan cir- 
cles show the stars with high resolution spectroscopy data from 
iMucciarelli et al.l ((2008i) (with numbers assigned as in this pa- 
per). The green cir cle shows the aperture with 20" radius that 
Pess ev etal] J2006) used for photometr y. Finally, the yellow cir- 
cle denotes the photometric aperture of iMucciarelli et al] (2006) 
with 90" radius. 



eluded that the first three stars outside of the SINFONI mosaic 
are cluster members and we included them in the final integrated 
spectrum for the cluster. 

Column 7 in Table [3] lists the stars, which we added to the 
final spectra for each globular cluster. 

4.3. Cluster light sampling 

Our study shows that with a reasonable amount of mosaic tiles 
and no more than 1.5 hours VLT observing time, we can sam- 
ple a significant fraction of the light from each of our sample 
of LMC GCs. We have seen that for well concentrated clusters 
there is no need for additional light sampling, other than a cen- 
tral mosaic, covering at least the half-light radius. The strategy 
of observing a mosaic on the cluster's centre and then a sequence 
of the brightest and closest stars within a certain radius, works 
very well in the case of sparse and not luminous clusters, as well 
as for those, which are rich, but not well concentrated. However, 
there are still some doubts about the cluster membership of the 
brightest stars around the intermediate age clusters, especially 
in the case of carbon rich stars. For this reason we would need 
high resolution spectroscopy to measure their radial velocities 
and chemical composition and compare with other RGB stars. 

Flux calibration of ground based spectroscopic data in the 
near-IR is particularly difficult, thus we cannot rely on direct 
luminosity estimates fr om our data. However, we can use the 
available photometry of Mucciarel li et al.1 d2006l) to estimate the 
approximate amount of sampled light in the intermediate age 
clusters. In this study the authors provide not only integrated 
magnitudes and colours, but also estimates of the M- and C-type 



stars contributions to the total cluster light in the /f-band. From 
there we know that the /f-band light of NGC 2162 is dominated 
by only one carbon rich star, which is responsible for ~ 60 % 
of the total cluster luminosity. This C-type star contributes with 
about 70 % to our final integrated cluster spectrum. From this 
we conclude that for this clus ter we are missing only ~ 10% of 
the luminosity, measured by Mu cciarelli et alJ (120061) . The sit- 
uation with the other two intermediate age clusters is similar. 
NGC 1806 is quite rich in stars in comparison with the other 
tw o clusters, as seen i n the K -band images in Fig.Q] According 
to IMucciarelli et all d2006l) 22 % of the total cluster light is 
coming from stars with K < 12. 3 m (AGB stars). About 77 % 
of it is due to four carbon stars. In our NGC 1806 integrated 
spectrum these four C-t ype stars account for ~6Q%. Following 
IMucciarelli et alJ d2006l) . 15 % of the light in NGC 2173 are due 
to one C-type star. In our spectrum this star is responsible for 
~20%. For the three old and metal poor clusters such detailed 
estimates of the contributions from different stellar phases are 
not available, thus we cannot make similar estimates for them as 
for the intermediate age clusters. Our central SINFONI mosaics 
cover more than one half-light radius. From this we conclude 
that the majority of the cluster light is covered. 

According to these rough estimates it is evident that we have 
reached our goal set up in Sect. Q] namely to obtain represen- 
tative luminosity weighted, integrated spectra in the /f-band for 
a group of clusters, having intermediate and old ages. The final 
spectra used for the scientific analysis in the following sections 
are shown in Fig. [6] The achieved signal-to-noise ratios are given 
in Table[5] The two clusters that have spectra dominated by car- 
bon rich giants, NGC 1 806 and NGC 2 1 62, visually seem to have 
a significantly lower S/N than the other clusters. However, as 
mentioned already in Sect. 13.51 this does not reflect the reality: 
these spectra are dominated by numerous absorption lines from 
carbon-based molecules, which the method for computing S/N 
interprets as noise. To properly estimate the S/N one would need 
synthetic spectra, including the full spectral synthesis. This is 
clearly not possible, as we are aiming to provide first templates 
that could allow computation of such spectra. We estimated the 
S/N for these two clusters based on the amount of integrated flux 
as compared to the other clusters in our sample and assigned to 
them a conservative lower limit of S/N > 80. 

5. Line strength indices in the A-band 

The stellar population modelling technique allows us to obtain 
detailed estimates of the properties of integrated stellar pop- 
ulations by measuring, for example, the strengths of selected 
absorption or emission features in the spectra and comparing 
the measured values to the predicted ones. This approach has 
shown to be very effective in the optical wavelength range 
and a well establi s hed system, the so called Lick/IDS sy stem 
dFaberet all [1981 IWorthev etal] fl99lt iTrager et alJ [19981) . is 
widely used. For the study of /T-band sp ectral features several in- 
dex definitions hav e been adopt e d (e.g . iFor ster Schreiber 2000; 
jFrogel et al] 120011: ISilva et al] 120081: iMarmol-Oueralto et al.l 
120081) . To measure the line stren gths of Na I and Ca I (see Fig. [6]) 
we used the index definitions of lFrogel et al] d200l[). and for the 
streng th of 12 CO (2-0) the d efinitions of Marmol-Oueralto et al.l 
(2008) and lMarastorJd2005h . 

The principle of measuring a near-IR index is the same as in 
the optical Lick system. The value is computed as the ratio of 
the flux in a central passband to the flux at the continuum level, 
measured in two pseudo-continuum passbands on both sides of 
the central passband. Due to the lack of a well defined contin- 
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Fig. 6. Final integrated spectra of the six LMC globular clusters, 
used for analysis. The spectra of NGC 1806 and NGC 2162 are 
clearly dominated by carbon rich stars, evident from the numer- 
ous carbon based absorption features (e.g. C2, CH, CN) and the 
flatter continuum shape. 



uum on the red side of the 12 CO (2-0) feature, the Deo index 
that lMarmol-Oueralto et al.l (|2008) defined, uses two continuum 
passbands on the blue side. This index measure the ratio be- 
tween the average fluxe s in the continuum and the absorption 
bands. Maraston (2005) uses a CO index, which measures the 
ratio of the flux densities at 2.37 and 2.22 /mi, based on the 
HST/NICMOS filters F237M and F222M. The index is com- 
puted in units of magnitudes and is normalised to Vega. This 
index reflects the strength not only of 12 CO (2-0), but also of 
the other CO absorption features in the range 2.3 - 2.4 pm (see 
Fig.©. 

Prior to index measurements we broadened our spectra to 
a spectral resolution of 6.9 A (FWHM) to match the resolu- 
tion of similar earlier studies of the K-bwd lig ht of early type 
galaxies and stars obtained w ith VLT/ISAAC dSilvaet alj |2008; 
Mar mol-Oueralto et alj2009). We measured the recession veloc- 
ities of the LMC globular clusters and stars with the IRAF task 
fxcor and corrected for them. We did not apply velocity dis- 
persion corrections to the indices due to the very low velocity 
dispersions of the clusters (< 10 km s , see previous section). 
The final index values are listed in Table[5] 

The current stellar populatio n models, wh ich include the 
near-IR wavelength range (e.g. Maraston 2005), have too low 
spectral resolution at these wavelengths (~ 200 A FWHM) to be 
able to give predictions for the weak and narrow features Na I 
and Ca I. Until more detailed models become available we can 
only empirically explore their dependance on the parameters of 
the globular clusters, derived from resolved light studies. 

In Fig. [7] we show the dependance of the Nal and Cal in- 
dex on the age of the clusters. We see that the Nal index in- 
creases with decreasing age of the cluster an d increasing me tal - 
licity. This result was already suggested by ISilva et al.l d2008). 
who find that the centres of early type galaxies in the Fornax 
cluster with signatures of a recent star formation, i.e. stronger 
HjS, have also stronger Nal indices. The Cal index seems to 
show similar behaviour as Dqo (discussed further in the text) 
as function of age, albeit with larger error bars. 

The 12 CO (2-0) absorption feature at 2.29 pm, which we de- 
scribed with the Deo index, is much stronger and broader. Stellar 
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Fig. 7. Near-IR index measurements plotted vs. the age of the 
clusters. Top panel: Nal, Bottom panel: Cal index. 

population model predictions for its strength exist and will be 
discussed in the following section in more detail. 



6. Comparisons with stellar population models 

One of the goals of this project was to verify the predictions of 
the current stellar population models in the n ear-IR wavelengt h 
range. Such models are the ones presented in Maraston (2005). 
The author provides a full set of SEDs (spectral energy distribu- 
tions) for different stellar populations with ages from 10 3 yr to 
15 Gyr and covering a range of metallicities. These SEDs extend 
up to 2.5 pm, but their spectral resolution in the near-IR wave- 
lengths is rather low, with one pixel covering 100 A. 

The model predictions include the CO index, and for 
completeness we measured the Deo index as defined by 
Mar mol-Oueralto et al.l d2008l) . In order to measure it for the 
models we had to interpolate the model SEDs linearly to a 
smaller wavelength step of 14 A. We have chosen this value 
after a few tests to check t he stability of the computation. 
Marmol- Oueralto et al.l d2008l) show that their Deo index defini- 
tion is very little dependent on instrumental or internal velocity 
broadening of the spectrum. However, the most extreme resolu- 
tion they have tested is ~ 70 A {FWHM), while the resolution 
of the model SEDs is ~ 200 A {FWHM). We have measured the 
Deo index in 10 stars from our sample with their nominal resolu- 
tion and when broadened to match the resolution of the models. 
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Name 


Mo T r Al 


r\. t r ai 


L>CO 


CU |magj 




(i) 


(2) 


(3) 


(4) 


(5) 


(6) 


NGC 1754 


0.14±0.33 


-0.11+0.74 


1.082+0.005 


-0.035+0.002 


55 


NGC 2005 


0.34±0.16 


0.18±0.35 


1.086+0.003 


0.017±0.001 


107 


NGC 2019 


0.32±0.25 


-0.45±0.51 


1.068+0.003 


-0.003+0.001 


103 


NGC 1806 


3.41±0.32 


1.27+0.65 


1.129±0.005 


0.026+0.001 


>80 


NGC 2162 


3.37±0.66 


0.67+1.35 


1.108+0.010 


0.013±0.002 


>80 


NGC 2173 


2.76±0.29 


2.19±0.63 


1.186+0.005 


0.123+0.002 


75 



Notes: (1) Cluster name, (2) Nal index, (3) Cal index, (4) Deo index, (5) CO index, (6) signal-to-noise ratio of the integrated spectrum. 



The broadened spectra have a Deo index value, which is weaker 
by 0.04 with respect to the unbroadened spectra. An offset with 
this size is not going to influence our general conclusions about 
the integrated spectra of LMC GCs. This is further supported by 
the CO index values of the GCs, which exhibit the same rela- 
tive values compared to the models. The CO index covers a very 
large wavelength range and is little dependent on the spectral res- 
olution. Thus we decided to measure the Deo index at the nomi- 
nal resolution of our spectra (6.9 A FWHM) and treat the model 
predictions with caution. Once higher resolution models become 
available, this comparison should be repeated in a more quantita- 
tive way. In the following subsection we discuss the comparison 
between the models and the data. For clarity we divided the glob- 
ular clusters in two groups. In Sect. |6.1| we explore the three old 
and metal poor clusters NGC 1754, NGC 2005 and NGC 2019. 
In Sect. 16.21 we discuss the intermediate age globular clusters 
NGC 1806, NGC 2162 and NGC 2173. 



6.1. Old clusters 

Literature data of the integrated near-IR colours of our sample 
of old and metal poor clusters are shown in Fig. [8] with open 
circles. The large open symbols represent the colours, derived 
from 100" radius aperture photometry, the small open symbol s 
from 20" radius aperture photometry from Pess ev etldl (|2006). 
With blue lines w e have overplotted SSP model predictions from 
Maraston (2005]) for a Kroupa IMF and blue horizontal branch 
morphology. The metallicity which is closest to the one, derived 
for our GCs sample is denoted with a solid blue line, [Z/H] = - 
1.35. In the top panel the data agree reasonably well with the 
model (J - K) colour. In the bottom panel the (H - K) colours 
show a larger spread. 

In Fig. [9] we show the comparison between model predic- 
tions and ind ex measurements for 12 CO (2-0). For the top panel 
we used the iMarastonl (|2005j) CO inde x, while for the bottom 
panel we used the D co index defined in Marmol-Oueralto et al. 
(2008). For stellar populations with an age of more than 3 Gyr 
the near-IR A"-band light is dominated by RGB stars, whose con- 
tribution stays approximately constant over large time scales. 
This is reflected in the stellar population models in Fig.|£l where 
the CO index remains almost constant at a given metallicity for 
ages > 3 Gyr. Our observational data fit reasonably well with the 
model predictions, despite the slightly lower index values in the 
LMC GCs with respect to the [Z/H] =-1.35 line. A possible rea- 
son for this, including the bluer colours of the globular clusters 
with respect to the models, may arise from the fact that the mod- 
els use [Z/H] to describe the metallicity, which includes not only 
iron but also other heavy elements. The literature data, which we 
used, estimate the metallicity only based on [Fe/H]. If the glob- 
ular clusters in our sample follow similar chemical trends as the 



old LMC globular clusters discussed in lMucciarelli et ail ((2009) 
we might expect a better agreement. 

6.2. Intermediate age clusters 

We divided the intermediate age GCs sample into two sub- 
samples, according to the SWB type of the clusters. NGC 1806 
and NGC 2162 have SWB t ype V and thus t heir age is esti- 
mated to be around 1.1 Gyr dFrogel et al.lll990l) . NGC 2173 has 
SWB between V and VI and an age of approximately 2 Gyr 
(Frog el et alj fl990). Based on colour-magnitude diagram meth- 
ods different authors give slightly higher or lower ages for these 
clusters. In order to be co nsistent wi t h the a ge calibration of the 
SSP models presented by Maraston (2005) we have decided to 
use the ages based on SWB types, as was done in these models. 

In Fig. [8] a compilation of literature photometric data for 
the intermediate age GCs in our sample are shown. According 
to their mean metallicity the closest model is shown with 
the dotted lines and has [Z/H] = -0.33. The integrated 
near-IR colours within an aperture with radius of 90" from 
iMucciarelli et a l. (2006) are shown with solid circles. The pho- 
tometry of Pe ssev etldl (120061) with different aperture sizes is 
shown with open circles. 

Comparing the large filled and open symbols in Fig. [8] we 
see evidence for disagreement between t he two studies, de- 
spite the fact that IMucciarelli et alj d2006l) do not provide er- 
ror bars. We also se e that according to the different apertures of 
iPessev et al.1 d2006) the intermediate age clusters exhibit colour 
gradients of about 0.2™ in (J - K) and (H - K). However, not 
always in the same direction. In (J - K) the clusters become 
systematically bluer with increasing radius, while in (H - K) 
NGC 2162 gets redder. In general colour gradients in globular 
clusters can be explained by mass segregation, which makes the 
most massive, and thus the most evolved, stars to be concentrated 
towards the centre. However, this effect i s not expected to be that 
large, in comparison with the values that Pessev et al. (2006) re- 
port. The A"-band magnitudes of IPessev et al.l (|2006) are fainter 
by appro ximately one magnitude as compared to the ones mea- 
sured bv lMucciarelli et all I d2006l) for NGC 1806 and NGC 2162 
for a similar aperture and the (J - K) colours are bluer. These 
effects could be either due to an ov erestimated LMC field de- 
contamination, if IPessev et al.1 (120061) removed the reddest stars, 
or und erestimation of the field in the case of IMucciarelli et ail 
(120061) . We can check this hypothesis using NGC 2162, where 
the A"-band light of the globular cluster is dominated by one 
single carbon rich star. By adding its A"-band magnitude, taken 
from the 2MASS catalogue, to the integrated (over 10 0") mag- 
nitude of the cluster, taken from IPessev et al.l d2006l) . we ob- 
tain a final AT-band ma gnitude in much better agreement with 
IMucciarelli et al.l ((2006). Field carbon stars are not observed fre- 
quently at the location of this cluster in the LMC, thus there is a 
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Fig. 8. Comparison of clusters' integrat ed colours with SSP 
model predictions from Maraston] d2005l) . The filled coloured 
circles (colour coding as listed in the upper right corner of the top 
plot) correspond to the (J-K) and ( H — K) colours derived form 
the 90" radius aperture photometry of Mu cciarelli et alj d2 006). 
The c oloured open circles give the values, which iPessev et alj 
(2006) derived from two different size apertures: small circles 
for r = 20", large circles for r = 100". The blue thick lines with 
different styles give the model predictions as a function of the 
metallicity and the age. 



Fig. 9. A comparison of model predictions of Maraston (2005) 
with our LMC GCs data (coloured filled circles). Model metal- 
licities are given with dif ferent line styles Top panel: CO index, 
used by Maraston (200 51 to compute CO line s trengt h. Bottom 
panel: D C o, defined bv lMarmol-Oueralto et all d2008l) . On both 
plots the arrows show the changes in CO after adding a different 
amount of carbon star contribution to the NGC2173 spectrum. 
From top to bottom the carbon star fraction is 40 %, 50 %, 60 % 
and 70 % respectively. The red arrow shows the CO index when 
the spectrum of NGC 2173 is composed of 6 % of the spectrum 
in the 3 rd bin of lLancon & Mouhcind d2002l) . 



high probability that this star is a cluster me mber. Therefore the 
large colour gradients present in the work of Pes sev et"ail ([2006) 
and their fainter /T-band magnitudes are most probably due to 
their oversubtraction of the LMC field star contribution for these 
two clusters. 

Our comparison of the SSP model predictions with the 
observed CO indices of LMC globular clusters is shown in 
Fig. [9] For the intermediate age clusters on both panels the 
closest model metallicity to the data is denoted with the dot- 
ted line, [Z/H]= -0.33. NGC 2173 (age ~2 Gyr) is marked with 
a red filled circle and agrees reasonably well with the model 
predictions. However, the other two intermediate age clusters, 
NGC 1806 (light green filled circle) and NGC 2 162 (orange 
filled circle) have CO index values much lower than the model 
predictions, independent from the CO index definition. 

6.2.1 . Cluster light sampling and stochastic effects 

We have tested a few scenarios, which can provide possible 
explanations of the observed trends. The first scenario, which 



we have investigated is whether we are sampling well the to- 
tal cluster light. We consider this is not a problem following 
the estimations about the sampled cluster light during our ob- 
servations, made in Sect. I4.3I The second question is, are the 
two clusters stochastical ly sampling the IMF well? H ere the an- 
swer is "most likely". Lancon & Mouhcine (2000) show that 
the CO index is significantly influenced by stochastic fluctua- 
tions for intermediate age stellar populations with 10 3 Mq, but 
when the mass reaches 1O 4 M0 and above the fluctuations are 
much less prominent. Our intermediate age globular clusters 
have masses closer to the ~ 10 4 M (T) range (e.g. NGC 2162, 
iMcLaughlin & van der M arel 2005) and thus are not expected 
to exhibit large variations in the CO index. By building a su- 
per cluster from NGC 1806 and NGC 2162 we are reducing the 
stochastic effects further, and we still find the same behaviour of 
decreasing Deo index value. However, this finding needs to be 
further verified with larger samples of intermediate age clusters. 
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Fig. 10. (J - K) vs. Deo index of the globular clusters and stars 
in our sample. The squares stand for the central mosaics of the 
GCs, the diamonds show the additional bright stars in our sam- 
ple. The stars denoted with circles are carbon-rich stars. 



6.2.2. Different Carbon stars in LMC and Milky Way? 

The careful investigation of the integrated spectra of LMC in- 
termediate age GCs showed us that the presence of C-type stars 
does not only influence the overall colours of the clusters, as 
discussed above, but also their CO index. Fig. [10] where we 
plotted the Dco index as a function of the ( J - K) colour for 
all of the objects in our observational sample, shows a general 
trend of increasing 12 CO (2-0) index with the colour becoming 
redder. This is typical for oxygen rich (or M-type) stars (e.g. 
Frogel et al. 1990). However, stars with {J-K) > 1.3 show a de- 
creasing 12 CO (2-0) line strength. These are exactly the carbon 
rich stars, separated in Fig. [2] with a slanted line, and denoted 
in Fig. [TOl with circles. Similar trends f or carbon stars of LMC 
globular clusters have been observed bv lFrogel et al.l (fl990). 

According to lFrogel et al.l (1 1990b the clusters, which harbour 
the brightest carbon stars are of SWB type V and VI. For these 
clusters, carbon stars contribute about 40 % of the total bolomet- 
ric light and this observation i s taken into coun t in the calibra- 
tions of the computed SEDs of iMarastonl J2005) models. Using 
the spectrum of NGC2173, which is very convenient for this 
kind of tests, because it is the one with a minimal carbon star 
contribution in our sample, we checked what is the effect of dif- 
ferent ratios of M-type to C-type stars on the CO line strength 
of the final integrated spectrum. The results are shown in Fig. [9] 
with black arrows, which start from the initial position of NGC 
2173 (red filled circle) and end at the age of clusters with SWB 
type V, i.e. 1.1 Gyr. The reason for this age scaling is that a larger 
carbon star contribution would mimic the spectrum of a younger 
cluster. The arrows show, from top to bottom, an increasing frac- 
tion of carbon star contribution to the final cluster light - 40 %, 
50 %, 60 % and 70 %, respectively originate from the C-type star . 
The same is valid for both the CO index used by Maraston ( 2005 ) 
and the Deo index defined by iMarmol-Oueralto et al.l Q2008). 
This simple experiment shows us that we are able to reproduce 
the lower observed CO index values in the integrated spectra of 
LMC globular clusters having SWB type V with an increasing 
fraction of the carbon star contribution. The same kind of tests, 
but performed with integrated colours, make the cluster redder, 



as expected when having an increased fraction of carbon rich 
stars. 

The models of Maraston (2005) include a careful treatment 
of the TP-AGB stellar phase, which is of great importance for 
stellar populations with ages between 0.3 and 2 Gyr, due to its 
very high luminosity. The empirical photometric calibration of 
the models has been done with the near-I R photometric data of 
LM C globular c l usters and AGB stars of iPersson et al.l (Il983l) 
and lFrogel et al.l dl990l) . Spectra of carbon -rich stars are also in- 
cluded . They come from the database of Lancon & Mou hcinel 
(120021) . which contains averaged C-type star spectra. These au- 
thors have obtained 21 spectra of carbon-rich luminous pulsat- 
ing variable stars in the Milky Way. However, due to the small 
temperature scales of the sample, they do not consider as jus- 
tified to have more than a few averaged bins. The carbon stars 
have been grouped according to their temperature, defined by 
their (H - K) colour. The first three bins contain averages of 6 
C-type star spectra each. Bins 4 and 5 contain the spectra of a 
single very red star, R Lep, near maximum and minimum light. 
The temperature of the stars is decreasing with increasing bin 
number. 

If carbon rich stars in the Galaxy and the LMC have differ- 
ent CO absorption strengths, then this would have a profound 
impact on the model predict ions. Differences in the CO index 
have been observed by e.g. iFrogel et al.l (1 19901) using narrow- 
band filters, where they compare the CO value of LMC clusters 
C-type stars and their counterparts in the Milky Way. The LMC 
stars have systematically weaker CO indices at a given colour. 
In Fig.[TT]we show the CO index values as a function of (J - K) 
colour, measured in C-type stars belonging to our sample (filled 
symbols), compa red to the CO index that we measured on the av- 
eraged spectra of Lancon & Mouhcine (2002). The least-squares 
linear fit to the LMC carbon rich stars is shown as a solid line: 



CO = 0.676 + 0.026 - 0.386 + 0.015 (J - K) 
D co = 1-668 + 0.103 - 0.321 + 0.058 (J - K) 



(2) 
(3) 



The Milky Way averaged spectra are marked with open sym- 
bols, with a bin number assigned to each. The dashed line is 
showing the linear least-squares fit to the averaged spectra: 



CO = 0.373 + 0.002 - 0.146 + 0.001 (J - K) 
D co = 1-268 + 0.007 - 0.042 + 0.003 (J - K) 



(4) 
(5) 



From this figure we see that the trends for C-type stars in 
the Milky Way and the LMC generally disagree. Moreover, stars 
with (J-K) > 1.6 have increasingly different CO index val- 
ues, which explains why our LMC GCs sample fits the model 
predictions for integrated colours, but not for the CO index. 

We made another test of including carbon star light to the 
spect rum of NGC 2173, but this time using the spectrum of bin 
3 of lLancon & Mouhcinel (120021) . The result is shown in Fig. [9] 
with a red arrow. For this case specifically the ratio of carbon 
star to the original spectrum was 6:4. Different ratios, as in the 
previous test, gave similar values. We see that the resultant CO 
index value increases by adding Milky Way carbon stars and fol- 
lows the trends predicted by the models, while the inclusion of 
LMC cluster carbon stars leads to a decrease of the index. We 
find this as an attractive explanation for the discrepancy between 
the models and the data. 

At this point the question arises why C-type stars in the LMC 
and the Galaxy have different CO line strengths at a given near- 
IR colour? We have to treat this issue with care, since not all of 
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Fig. 11. (J - K) vs. CO index in LMC and Milky Way carbon- 
rich stars. LMC stars are marked with filled circles. The lin- 
ear least-squares fit to the data is shown with a solid line 
(Eq.0. T he values fo r the M ilky Way averaged star spectra from 
Lancon & Mouhcin^ d2002l) are shown with open circles. The 
number next to each data point corresponds to the bin number. 
Their linear least-squares fit is shown with a dashed line (Eq.|5]l. 

the carbon stars in the Lancon & Mouhcine] d2002l) library dis- 
agree with the relation for the LMC C-stars. The averaged spec- 
trum from bin 2 agrees with the LMC relation. In the bottom 
panel of Fig. QT|we indicated with an asterisk another C-type 
star from the Milky Way, located in the Galac tic open cluster 
NGC 2477 ([Fe/H] = -0.02, (J-K) = 1.69 dHoudashelt et alj 
19921). A A"-band spectrum of this star has been published in 
Silva et al] d2008l) . however the authors did not discuss its prop- 
erties. Its colour and Deo index value are consistent with the 
relation for LMC C-type stars, given in Equation|2] 

AGB stars in populations with lower metallicity are more 
likely to become C-type. A star becomes C-type when the ra- 
tio of carbon to oxygen atoms in its atmosphere becomes larger 
than one. Then all the oxygen is locked up in CO molecules, 
and there is still some extra carbon in the atmosphere. For the 
metal-poor stars there is less oxygen in the atmosphere in the 
first place. Therefore these stars need to dredge up less carbon 
in order to overcome, in relative amount, the quantity of oxygen 
and so become C-type. The remaining carbon is used to form 
other molecules, like C2, CN and CH. In this way the CO index 
as we measure it decreases. Note that the overall line strength 
of the 12 CO(2-0) may not decrease significantly. For exam- 



ple see Fig.|6l where the 12 CO(2-0) features in NGC 1806 and 
NGC 2162 do not seem to be much weaker, than in NGC 2173. 
But the CO indices, which we use to describe the behaviour 
of this feature, take into account the continuum shape as well, 
which in the case of carbon stars is severely affected by the pres- 
ence of features like C2, CN or CH. In more metal rich stars the 
amount of oxygen atoms is higher and all of the available carbon 
is used to form CO molecules, which makes a stronger CO index 
in oxygen-rich stars. The different CO indices in Milky Way and 
LMC carbon-rich stars are due to a real change in the depth of 
the CO absorption features. 

Th e observed difference between the lLancon & Mou hcinel 
d2002l) library carbon stars in the Milky Way and the ones in 
our sample in the LMC has important implications on the stel- 
lar populations models. The photometric ca libration of the i nter- 
mediate age populations in the models of lMarastonl d2005l) has 
been performed with globular clusters in the LMC. However, for 
the spectral calibration, spectra of carbon-rich stars in the Milky 
Way have been used, which leads to inconsistent results. There 
might be a number of reasons for the observed inconsistencies, 
like different metallicities of the carbon stars or the phase of the 
pulsation in which they were observed. Our results clearly show 
the necessity of a better understanding of the properties of car- 
bon rich stars and the need of larger empirical stellar libraries to 
improve model predictions. 



7. Conclusions 

The goal of this project is to provide an empirical spectral li- 
brary in the near-IR for integrated stellar populations with ages 
> 1 Gyr, which will be used to test the current and calibrate 
the future stellar population models. In this paper we have pre- 
sented the first results from a pilot study of the /T-band spectro- 
scopic properties of a sample of six globular clusters in the LMC. 
To validate the observational strategy, data reduction and analy- 
sis me thods, we have selected from the catalogue of iBica et al.l 
d 19991) three out of 38 GCs with SWB type VII to represent the 
old (>10 Gyr) and metal poor ([Fe/H]~-1.4) population of the 
LMC, and three out of 71 clusters with SWB types V and VI to 
explore the properties of the intermediate age (1-3 Gyr) and 
more metal rich ([Fe/H] — 0.4) component of the population. 
For each cluster our integrated spectroscopy covers the central 
24" x 24" and in most of the cases we have sampled about half 
the light. However, in order to better sample bright AGB stars, 
which are the most important contributors to the integrated clus- 
ter light in the near-IR, we have observed up to 9 of the brightest 
stars outside the central mosaics, but still within the tidal radii 
of the clusters, that have near-IR colours and magnitudes consis- 
tent with bright red giants in the observed clusters. We obtained 
integrated luminosity weighted spectra for the six clusters, mea- 
sured the line strengths of Nal, Cal and 12 CO(2-0) absorption 
features in the A'-band and compare d the strength of 1 2 CO (2-0) 
with the stellar population models of lMarastonl d2005l) . 

The observing strategy to cover at least the central half-light 
radius with a number of SINFONI pointings showed to be an 
efficient way in sampling the near-IR light of old (> 10 Gyr) 
clusters. For the intermediate age and sparse clusters, which are 
dominated by just a few very bright stars, observing a central 
mosaic plus a number of the brightest stars in the vicinity of 
the cluster is a better choice for optimal cluster light sampling. 
The availability of high spectral resolution spectroscopy greatly 
helps the differentiation of the cluster member stars from the 
LMC field population. 
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In intermediate age clusters the largest amount of light orig- 
inates from oxygen (M-type) and carbon-rich (C-type) AGB 
stars. Different ratios of the contributions by these two types 
of stars can lead to significant changes in the near-IR 12 CO (2- 
0) line strength. According to our observations, when the C- 
type stars contribution peaks (at ~1 Gyr) the observed CO line 
strength is weak and then increases fast to reach its maximum 
for clusters with age ~2 Gyr. It is important to note that a weak 
line strength of 12 CO (2-0) does not mean that there is less CO in 
these clusters/stars. The indices, which are used to describe the 
line strengths, take also into account the continuum shape, which 
in the case of carbon-rich stars is severely affected by typical for 
this type of stars absorption features and thus the resultant index 
value is low. 

T he comparison o f our data with the stellar population mod- 
els of Maraston ( 2005), in terms of CO line strength, shows a dis- 
agreement for the youngest clusters in our sample. For clusters 
with age ~ 1 Gyr the models predict maximal CO line strength, 
while we observe the opposite: the CO strength is significantly 
weaker. At the same time, literature data of the integrated colours 
of the clusters are consistent with these models. We explain these 
discrepancies as due to the different origin of the C-type stars 
used to calibrate the models and the ones in our data sample. The 
stars, used for model calibration, are Milky Way carbon stars, 
while our carbon stars are born in the LMC. We support this sce- 
nario with Fig.QT| where we show that carbon rich stars in the 
Milky Way and LMC, which have similar (J - K) colour, have 
very different CO line strengths. This deserves further investi- 
gatio n and hopefully the next generation of carbon star models 
(e.g. lAringer et al.l l2009) will help us to elucidate whether a sys- 
tematic effect of the metallicity on CO indices is expected, or 
the found discrepancy is due to the small sample of individual 
observations of variable stars in existing libraries. 

The near-IR Na I index shows a dependance on the age of the 
clusters - it is decreasing with an increasing age. The combina- 
tion between optical and near-IR spectral indices seems to offer 
possibilities to break the age-metallicity degeneracy, but more 
accurate and detailed stellar population models are necessary in 
the near-IR wavelength range. 

These models are of paramount importance, when study- 
ing the spatially resolved stellar populations of nearby galaxies. 
Adaptive optics assisted observations allow for the best correc- 
tion of the Earth's atmosphere perturbing effects when observing 
in the near-IR. First attempts to explore g alaxy evolution via spa- 
tially resolved near-IR spectroscopy (e.g . lLvubenova et alT2 008; 
iDavidge et al.l2008l : lNowak et alj2 008) have shown the need for 
a better understanding of the properties of stellar populations in 
this wavelength range. The availability of detailed and reliable 
stellar population models in the near-IR will open up a new win- 
dow to the exploration of galaxy formation and evolution. 

The final integrated spectra of the six LMC globular clus- 
ters will be made available via the Strasbourg astronomical Data 
Center (CDS). 
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